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SYSTEM OF SUPER SUPER DECOUPLED LOADFLOW CALCULATION FOR 

ELECTRICAL POWER SYSTEM 

FIELD OF THE INVENTION 

[001] The present invention relates to a method of load flow calculation in power flow control 
and voltage control for an electrical power system. 

BACKGROUND OF THE INVENTION 

[002] The present invention relates to power-flow/voltage control in utility/industrial power 
networks of the types including many power plants/generators interconnected through 
transmission/distribution lines to other loads and motors. Each of these components of the power 
network is protected against unhealthy or alternatively faulty, over/under voltage, and/or over 
loaded damaging operating conditions. Such a protection is automatic and operates without the 
consent of power network operator, and takes an unhealthy component out of service by 
disconnecting it from the network. The time domain of operation of die protection is of the order 
of milliseconds. 

[003] The purpose of a utility/industrial power network is to meet the electricity demands of its 
various consumers 24-hours a day, 7-days a week while maintaining the quality of electricity 
supply. The quality of electricity supply means the consumer demands be met at specified voltage 
and frequency levels without over loaded, under/over voltage operation of any of the power 
network components. The operation of a power network is different at different times due to 
changing consumer demands and development of any faulty/contingency situation. In other words 
healthy operating power network is constantly subjected to small and large disturbances. These 
disturbances could be consumer/operator initiated, or initiated by overload and vmder/over voltage 
alleviating functions collectively referred to as security control functions and various optimization 
functions such as economic operation and minimization of losses, or caused by a 
fault/contingency incident. 

[004] For example, a power network is operating healthy and meeting quality electricity needs 
of its consumers. A fault occurs on a line or a transformer or a generator which faulty component 
gets isolated from the rest of the healthy network by virtue of the automatic operation pf its 
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protection. Such a disturbance would cause a change in the pattern of power flows in the network, 
which can cause over loading of one or more of the other components and/or over/under voltage 
at one or more nodes in the rest of the network. This in turn can isolate one or more other 
components out of service by virtue of the operation of associated protection, which disturbance 
can trigger chain reaction disintegrating the power network. 

[005] Therefore, the most basic and integral part of all other functions including optimizations 
in power network operation and control is security control. Security control means controlling 
power flows so that no component of the network is over loaded and controlling voltages such 
that there is no over voltage or under voltage at any of the nodes in the network following a 
disturbance small or large. As is well knovra, controlling electric power flows include both 
controlling real power flows which is given in MWs, and controlling reactive power flows which 
is given in MVARs. Security control functions or alternatively overloads alleviation and 
over/under voltage alleviation functions can be realized through one or combination of more 
controls in the network. These involve control of power flow over tie line connecting other utility 
network, turbine steam/water/gas input control to control real power generated by each generator, 
load shedding function curtails load demands of consumers, excitation controls reactive power 
generated by individual generator which essentially controls generator terminal voltage, 
transformer taps control coimected node voltage, switching in/out in capacitor/reactor banks 
controls reactive power at the connected node. 

^^1^6^ Control of an electrical power system involving power-flow control and voltage control 
commonly is performed according to a process shown in fig. 5. The various steps entail the 
following. 

Step- 10: Obtain on-line/simulated readings of open/close status of all switches and circuit 
breakers, and read data of maximum and minimum reactive power generation capability 
limits of PV-node generators, maximum and minimum tap positions limits of tap 
changing transformers, and maximum power carrying capability limits of transmission 
lines, transformers in the power network, or alternatively read data of operating limits of 
power network components. 
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Step-20: Obtain on-line readings of real and reactive power 
assignments/schedules/specifications/settings at PQ-nodes, real power and voltage 
magnitude assignments/schedules/specifications/settings at PV-nodes and transformer 
turns ratios. These assigned/set values are controllable and are called controlled 
variables/parameters . 

Step-30: Resulting voltage magnitudes and angles at power network nodes, power flows through 

various power network components, reactive power generations by generators and turns 

ratios of transformers in the power network are determined by performance of loadflow 

calculation, for the assigned/set/given/known values from step-20 or previous process 

cycle step-60, of controlled variables/parameters. 
Step-40: The results of Loadflow calculation of step-30 are evaluated for any over loaded power 

network components like transmission lines and transformers, and over/under voltages 

at different nodes in the power system 
Step-50: If the system state is acceptable implying no over loaded transmission lines and 

transformers and no over/under voltages, the process branches to step-70, and if 

otherwise, then to step-60 
Step-60: Changes the controlled variables/parameters set in step-20 or as later set by the previous 

process cycle step-60 and returns to step-30 
Step-70: Actually implements the corrected controlled variables/parameters to obtain 

secure/correct/acceptable operation of power system 

[007] Overload and under/over voltage alleviation functions produce changes in controlled 
variables/parameters in step-60 of Fig.5. In other words controlled variables/parameters are 
assigned or changed to the new values in step-60. This correction in controlled 
variables/parameters could be even optimized in case of simulation of all possible imaginable 
disturbances including outage of a line and loss of generation for corrective action stored and 
made readily available for acting upon in case the simulated disturbance actually occurs in the 
power network. In fact simulation of all possible imaginable disturbances is the modem practice 
because corrective actions need be taken before the operation of individual protection of the 
power network components. 
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[008] It is obvious that loadflow calculation consequently is performed many times in real-time 
operation and control environment and, therefore, efficient and high-speed loadflow calculation is 
necessary to provide corrective control in the changing power system conditions including an 
outage or failure of any of the power network components. Moreover, the loadflow calculation 
must be highly reliable to yield converged solution under a wide range of system operating 
conditions and network parameters. Failure to yield converged loadflow solution creates blind 
spot as to what exactly could be happening in the network leading to potentially damaging 
operational and control decisions/actions in capital-intensive power utilities. 

[009] The power system control process shqwn in Fig. 5 is very general and elaborate. It 
includes control of power-flows through network components and voltage control at network 
nodes. However, the control of voltage magnitude at connected nodes within reactive povver 
generation capabilities of electrical machines including generators, synchronous motors, and 
capacitor/inductor banks, and within operating ranges of transformer taps is normally integral part 
of load flow calculation as described in "LTC Transformers and MVAR violations in the Fast 
Decoupled Load Flow, IEEE Trans., PAS-101, No.9, PP. 3328-3332, September 1982." If 
under/over voltage still exists in the results of load flow calculation, other control actions, manual 
or automatic, may be taken in step-60 in the above and in Fig.5. For example, under voltage can 
be alleviated by shedding some of the load connected. 

[010] The prior art and present invention are described using the following symbols and terms: 



Ypq Gpq+jBpq 

j^jL^ Ypp = Gpp+jBpp 



yp = gp+jbp 

Vp = ep+jfp=VpZ0, 
AGp, Avf~ 
Acp, Afp 

Pp+jQp 
APp+jAQp 

RPp+jRQp 

On 



(p-q) th element of nodal admittance matrix without shunts 
p-th diagonal element of nodal admittance matrix without shunts 



total shxmt admittance at any node-p 

CO 



mplex voltage of any node-p -j^ _ ^ ^j'^ c /d?^ 4 o^pi^c^ ^cAa-ic-W^/^ 



voltage angle, magnitude corrections 
real, imaginary components of voltage corrections 
net nodal injected power calculated 
nodal power residue or mismatch 
modified nodal power residue or mismatch 
rotation or transformation angle 
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[RP] 



[RQ] 



: vector of mnHififH rpt^i pn^f>r r^ffi^iii" f r m bmatch 
: vector of modified reactive power residue or mismatch 



[Ye] 

[YV] 



m 
k 

n=m+k+ 1 
q>p 

[ ] 



: gain matrix of the P-9 loadflow sub-problem defined by eqn. ( 1) 
: gain matrix of the Q-V lo adflow sub-problem defined by eqn. (2) 
: number of PQ-nodes 
: number of P V-nodes 
: total number of nodes 

: q is the node adjacent to node-p excluding the case of q=p 
: indicates enclosed variable symbol to be a vector or a matrix 



: Limiting Rotation Angle, -36° for prior art, -48° for invented models 
PQ-node: load-node, where, Real-Power-P and Reactive-Power-Q are specified 
PV-node: generator-node, where, Real-Power-P and Voltage-Magnitude-V are specified 
Ypq' = Gpq' + jBpq': rotated (p-q) th element of nodal admittance matrix without shunts 
Ypp' = Gpp' + jBpp': rotated p-th diagonal element of nodal admittance matrix without shunts 
APp'= APpCos<Dp + AQpSin(Dp: rotated or transformed real power mismatch 
AQp'= AQpCos0p - APpSinOp: rotated or transformed reactive power mismatch 
Loadflow Calculation: Each node in a power network is associated with four electrical quantities, 

which are voltage magnitude, voltage angle, real power, and reactive power. 
The loadflow calculation involves calculation/determination of two unknown 
electrical quantities for other two given/specified/scheduled/set/known 
electrical quantities for each node. In other words the loadflow calculation 
involves determination of unknown quantities in dependence on the 
given/specified/scheduled/ set/known electrical quantities, 
a set of equations describing the physical power network and its operation for 
the purpose of loadflow calculation. The term 'loadflow model' can be 
alternatively referred to as 'model of the power network for loadflow 
calculation'. The process of writing Mathematical equations that describe 
physical power network and its operation is called Mathematical Modeling. If 
the equations do not describe/represent the power network and its operation 
accurately the model is inaccurate, and the iterative loadflow calculation 
method could be slow and unreliable in yielding converged loadflow 



Loadflow Model 
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calculation. There could be variety of Loadflow Models depending on 
organization of set of equations describing the physical power network and its 
operation, including Decoupled Loadflow Models, Super Decoupled 
Loadflow Models, Fast Super Decoupled Loadflow (FSDL) Model, and Super 
Super Decoupled Loadflow (SSDL) Model. 
Loadflow Method: sequence of steps used to solve a set of equations describing the physical 

power network and its operation for the purpose of loadflow calculation is 
called Loadflow Method, which term can alternatively be referred to as 
^loadflow calculation method' or 'method of loadflow calculation'. One word 
for a set of equations describing the physical power network and its operation 
is: Model. In other words, sequence of steps used to solve a Loadflow Model 
is a Loadflow Method. The loadflow method involves definition/formation of 
a loadflow model and its solution. There could be variety of Loadflow 
Methods depending on a loadflow model and iterative scheme used to solve 
the model including Decoupled Loadflow Methods, Super Decoupled 
Loadflow Methods, Fast Super Decoupled Loadflow (FSDL) Method, and 
Super Super Decoupled Loadflow (SSDL) Method. All decoupled loadflow 
methods described in this application use either successive (19, IV) iteration 
^ ^ ^ , scheme or simultaneous (IV, 10) iteration scheme, defined in the following. 



include a class of methods known as decoupled loadflow. This class of methods consists of 
decouled loadflow and super decoupled loadflow methods including Fast Super Decoupled 
Loadflow method. However, functional forms of different elements of the prior art Super 
Decoupled Loadflow: Fast Super Decoupled Loadflow (FSDL) model defined by system of 
equations (1) and (2) will be given below before description of steps of the prior art loadflow 
calculation method. The prior art FSDL model is very sensitive, in terms of iterations required to 
achieve converged loadflow calculation, to the restricted value of the rotation angle appUed to 
complex power mismatch in terms of mismatch in real and reactive power flowing in through 
each of PQ-nodes. Moreover, the presence of twice the value of transformed network shunt in a 
diagonal element of the gain matrix [YV] for certain power network causes slow convergence 
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taking increased number of iterations to converge to a solution, and therefore required increased 
calculation time. 

[012] The aforesaid class of Decoupled Loadflow models involves a system of equations for the 
separate calculation of voltage angle and voltage magnitude corrections. Each decoupled model 
comprises a system of equations (1) and (2) differing in the definition of elements of [RP], [RQ], 
[Y9] and [YV], 



[RP] = [YG] [AG] 
[RQ] = [YV] [AV] 



(1) 

(2) 



[013] A decoupled loadflow calculation method involves solution of a decoupled loadflow 
model comprising system of equations (1) and (2) in an iterative manner. Commonly, successive 
(19, IV) iteration scheme is used for solving system of equations (1) and (2) aitemately with 
intermediate updating. Each iteration involves one calculation of [RP] and [AG] to update [G] and 
then one calculation of [RQ] and [AV] to update [V]. The sequence of equations (3) to (6) depicts 
the scheme. 



[AG] = [YG] ■' [RP] 
[G] = [G] + [AG] 
[AV] = [YV] [RQ] 
[V] = [V] + [AV] 



(3) 
(4) 

I 

(5) 
(6) 



[014] The elements of [RP] and [RQ] for PQ-nodes are given by equations (7) to (10). 

RPp = (APpCosOp + AQpSin<Dp)A^p = APp Wp (7) 

RQp = (-APpSinOp + AQpCosOp)A^p = AQp Wp (8) 

CosOp = Absolute (Bpp / v'(Gi7+B^)) ^ Cos (-36°) (9) 



SinOp = -Absolute (Gpp / V(G^^+B^) > Sin (-36°) (10) 

[015] A description of Super Decoupling principle and the prior art FSDL model is given in, 
"Fast Super Decoupled Loadflow", ffiE proceedings Part-C, Vol.139, No.l, pp. 13-20, January 
1992. 
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[016] Fast Super Decoupled Loadflow (FSDL) model consists of equations (3) to (16). 



RPp = APp / (KpVp) 



Ye 



pq 



-Y 



pq 



(Bpq + 0.9(Ypq-Bpq)) 

•B. 



'pq 



YVp,- 



-Y 



pq 



-for branch r/x ratio < 2.0 
-for branch r/x ratio > 2.0 

-for branches connected between two PV-nodes or 
a PV-node and the slack-node 

-for branch r/x ratio < 2.0 



-(Bpq + 0.9(Ypq-Bpq)) -for branch r/x ratio > 2.0 



YGpp = I- ye 



q>p 



pq 



and 



YVpp = .2bo' + Z-YV 



q>p 



pq 



(11) 



(12) 



(13) 
(14) 



bp' = bpCosOp 



or 



bp' — bp 



Kp = Absolute (Bpp/YOpp) 



(15) 
(16) 



[017] Branch admittance magnitude in (12) and (13) is of the same algebraic sign as its 
susceptance. Elements of the two gain matrices differ in that diagonal elements of [YV] 
additionally contain the b' values given by equation (15) and in respect of elements 
corresponding to branches connected between two PV-nodes or a PV-node and the slack-node. In 
two simple variations of the FSDL model, one is to make YVpq=Y9pq and the other is to make 

Y9pq=YVpq. 

[018] The steps of loadflow calculation method FSDL are shown in the flowchart of Fig. 1. 
Referring to the flowchart of Fig.l, different steps are elaborated in steps marked with similar 
letters in the following. The words "Read system data" in Step-a correspond to step- 1 0 and step- 
20 in Fig. 5, and step- 14, step-20, step-32, step-44, step-50 in Fig. 6. All other steps in the 
following correspond to step-30 in Fig.5, and step-60, step-62, and step-64 in Fig. 6. 

a. Read system data and assign an initial approximate solution. If better solution estimate is 
not available, set all the nodes voltage magnitudes and angles equal to those of the slack- 
node. This is referred to as the slack-start. 



b. Form nodal admittance matrix, and Initialize iteration count ITRP = ITRQ= r = 0. 
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Compute sine and cosine of rotation angles using equations (9) and (10), and store them. 
If CosOp < Cos (-36°), set CosOp = Cos (-36°) and SinOp = Sin (-36°). 
Form (m+k) x (m+k) size gain matrices [YG] and [YV] of (1) and (2) respectively each in 
a compact storage exploiting sparsity, using equations (12) to (15). In [YV] matrix, 
replace diagonal elements corresponding to PV-nodes by very large value, say, 10.0^^. 
Factorize [Y0] and [YV] using the same ordering of nodes regardless of node-types and 
store them using the same indexing and addressing information. 

Compute residues [AP] at PQ- and PV-nodes and [AQ] at only PQ-nodes. If all are less 
than the tolerance (e), proceed to step-m. Otherwise follow the next step. 
Compute the vector of modified residues [RP] using (7) for PQ-nodes, and using (1 1) and 
(16) for PV-nodes. 

Solve (3) for [AG] and update voltage angles using, [9] = [9] + [A9]. 

Set voltage magnitudes of PV-nodes equal to the specified values, and Increment the 

iteration count ITRP=1TRP+1 and r=(ITRP+ITRQ)/2. 

Compute residues [AP] at PQ- and PV-nodes and [AQ] at PQ-nodes only. If all are less 
than the tolerance (e), proceed to step-m. Otherwise follow the next step. 
Compute the vector of modified residues [RQ] using (8) for only PQ-nodes. 
Solve (5) for [AV] and update PQ-node magnitudes using [V] = [V] + [AV]. While 
solving equation (5), skip all the rows and columns corresponding to PV-nodes. 
Calculate reactive power generation at PV-nodes and tap positions of tap-changing 
transformers. If the maximum and minimum reactive power generation capability and 
transformer tap position limits are violated, implement the violated physical limits and 
adjust the loadflow solution by the method described in "LTC Transformers and MVAR 
violations in the Fast Decoupled Load Flow, IEEE Trans., PAS-101, No.9, PP. 3328-3332, 
September 1982". 

Increment the iteration count 1TRQ=1TRQ+1 and r=(ITRP+lTRQ)/2, and Proceed to 
step-d. 

From calculated values of voltage magnitude and voltage angle at PQ-nodes, voltage angle 
and reactive power generation at PV-nodes, and tap position of tap changing transformers, 
calculate power flows through power network components. 
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{^^019]) In super decoupled loadflow models [YG] and [YV] are real, sparse, symmetrical and built 
only from network elements. Since they are constant, they need to be factorized once only at the 
start of the solution. Equations (1) and (2) are to be solved repeatedly by forward and backward 
substitutions. [Y9] and [YV] are of the same dimensions (m+k) x (m+k) when only a 
row/column of the slack-node or reference-node is excluded and both are triangularized using the 
same ordering regardless of the node-types. For a row/column corresponding to a PV-node 
excluded in [YV], use a large diagonal to mask out the effects of the off-diagonal terms. When 
the PV-node is switched to the PQ-type, removing the large diagonal reactivates the row/coliunn 
corresponding to a switched PV-node to PQ-node type. This technique is especially useful in the 
treatment of PV-nodes in the gain matrix [YV]. 

(^J02^ The convergence of the prior art FSDL method is very sensitive to the value of the 
restriction applied to the rotation angle. The best possible convergence from non-linearity 
consideration could be achieved by restricting rotation angle to maximum of -36 . However, 
when large Resistance (R)/Reactance (X) ratio branch is present in the network without creating 
non-linearity problem, it takes large number of iterations to converge. Moreover, the presence of 
twice the transformed value of network shunts in the diagonal elements of the gain matrix [YV], 
causes it to take increased number of iterations in case of certain power networks. These 
problems are overcome by formulating power flow equations such that transformed values of 
known/given/specified/scheduled/set quantities appears in the diagonal elements of the gain 
matrix [YV] as described in the following. , 

f(fAM[J>^\j\^^ J ^ ("summary of the INVENfltmp $wp|)(?./v^ ('•^^''•^ " ^ ^ ^ 

[021] It is a primary object of the present invention to improve convergence and efficiency of 
the prior art Fast Super Decoupled Loadflow calculation method under wide range of system 
operating conditions and network parameters for use in power flow control and voltage control in 
the power system. A further object of the invention is to reduce computer storage/memory or 
calculating volume requirements. 

[022] The above and other objects are achieved, according to the present invention, with any 
one system of Super Super Decoupled Loadflow, X'X' -version (SSDL-X'X') and its many 
manifestations, calculation for Electrical Power System. In context of voltage control, the 
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inventive system of SSDL calculation for Electrical Power system consisting of plurality of 
electromechanical rotating machines, transformers and electrical loads connected in a network, 
each machine having a reactive power characteristic and an excitation element which is 
controllable for adjusting the reactive power generated or absorbed by the machine, and some of 
the transformers each having a tap changing element, which is controllable for adjusting turns 
ratio or alternatively terminal voltage of the transformer, said system comprising: 

means defining and solving loadflow model of the power network characterized by 
inventive SSDL model for providing an indication of the quantity of reactive 
power to be supplied by each generator including the slack/reference node 
generator, and for providing an indication of turns ratio of each tap-changing 
transformer in dependence on the obtained-online or given/specified/set/known 
controlled network variables/parameters, and physical limits of operation of the 
network components, 

machine control means connected to the said means defining and solving loadflow model 
and to the excitation elements of the rotating machines for controlling the 
operation of the excitation elements of machines to produce or absorb the amount 
of reactive power indicated by said means defining and solving loadflow model in 
dependence on the set of obtained-online or given/specified/set controlled network 
variables/parameters, and physical limits of excitation elements, 

transformer tap position control means connected to the said means defining and solving 
loadflow model and to the tap changing elements of the controllable transformers 
for controlling the operation of the tap changing elements to adjust the turns ratios 
of transformers indicated by the said means defining and solving loadflow model 
in dependence on the set of obtained-online or given/specified/set controlled 
network variables/parameters, and operating limits of the tap-changing elements. 

[023] The method and system of voltage control according to the preferred embodiment of the 
present invention provide voltage control for the nodes connected to PV-node generators and tap 
changing transformers for a network in which real power assignments have already been fixed. 
The said voltage control is realized by controlling reactive power generation and transformer tap 
positions. 
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[024] The inventive system of Super Super Decoupled Loadflow (SSDL) calculation can be 
used to solve a model of the Electrical Power System for voltage control. For this purpose real 
and reactive power assignments or settings at PQ-nodes, real power and voltage magnitude 
assignments or settings at PV-nodes and transformer turns ratios, open/close status of all circuit 
breaker, the reactive capability characteristic or curve for each machine, maximum and minimum 
tap positions limits of tap changing transformers, operating limits of all other network 
components, and the impedance or admittance of all lines are supplied. A decoupled loadflow 
system of equations (1) and (2) is solved by an iterative process until convergence. During this 
solution the quantities which can vary are the real and reactive power at the reference/slack node, 
the reactive power set points for each PV-node generator, the transformer transformation ratios, 
and voltages on all PQ-nodes nodes, all being held within the specified ranges. When the iterative 
process converges to a solution, indications of reactive power generation at PV-nodes and 
transformer turns-ratios or tap-settings are provided. Based on the known reactive power 
capability characteristics of each PV-node generator, the determined reactive power values are 
used to adjust the excitation current to each generator to establish the reactive power set points. 
The transformer taps are set in accordance with the turns ratio indication provided by the system 
of loadflow calculation. 

[025] For voltage control, system of SSDL calculation can be employed either on-line or off- 
line. In off-line operation, the user can simulate and experiment with various sets of operating 
conditions and determine reactive power generation and transformer tap settings requirements. A 
general-purpose computer can implement the entire system. For on-line operation, the load flow 
calculation system is provided with data identifying the current real and reactive power 
assignments and transformer transformation ratios, the present status of all switches and circuit 
breakers in the network and machine characteristic curves in steps- 10 and -20 in Fig. 5, and steps 
12, 20, 32, 44, and 50 in Fig 6 described below. Based on this information, a model of the system 
based on gain matrices of invented load flow calculation systems provide the values for the 
corresponding node voltages, reactive power set points for each machine and the transformation 
ratio and tap changer position for each transformer. 

[026J The present inventive system of loadflow calculation for Electrical Power System consists 
of, a Super Super Decoupled Loadflow: X'X'-version (SSDL-X'X') method characterized in that 
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1) modified real power mismatch at any PQ-node-p is calculated as RPp = [APp' + (GppV Bpp') 
^Qp'] / Vp^ which takes different form for different manifestation of the SSDL-X'X' method, 2) 
the presence of transformed values of known/given/specified/scheduled/set quantities in the 
diagonal elements of the gain matrix [YV] of the Q-V sub-problem, and 3) transformation angles 
are restricted to maximum of -48'', and these inventive loadflow calculation steps together yield 
some processing acceleration and consequent efficiency gains, and are each individually 
inventive, and 4) except in the manifestation SSDL-BX, modified real power mismatches at PV- 
nodes are determined as RPp = APp/(KpVp) in order to keep gain matrix [YG] symmetrical. If the 
value of factor Kp=l, the gain matrix [YB] becomes unsymmetrical in that elements in the rows 
corresponding to PV-nodes are defined without transformation or rotation applied, as Yepq= -Bpq. 

[027] An alternate inventive system of loadflow calculation for Electrical Power System 
consists of, a Super Super Decoupled Loadflow: BGX'-version (SSDL-BGX') method, using 
simultaneous (IV, 10) iteration scheme, characterized in that 1) the modified real power 
mismatch at any PV-node-p is calculated as: 

RPp = [(APpA^p) - 1 GpqAV J / Vp .... _ ( ?0 <^ 

and 2) the modified real power mismatch at any PQ-node-p is calculated as: 

RPp = [(APpA'p) - Z GpqA Vq - (gp' A Vp)] / Vp or ^ j) 

RPp= [(APpA^p).lGpqAVq]/Vp {-y^? 

q=l 

[028] A further alternate inventive system of loadflow calculation for Electrical Power System 
consists of, a Super Super Decoupled Loadflow:X'GpvX' -version (SSDL-X'GpvX') method, using 
simultaneous (IV, 10) iteration scheme, characterized in that 1) the modified real power 
mismatch at any PV-node-p is calculated as: 

RPp= [(APpA^p) - Z Gpq AV,] / (Kp* Vp) _ (^?^ 

and, 2) the modified real power mismatch at any PQ-node-p is calculated as: 



RP, = [{[APp' + (Gp,VBpp')AQp'] / Vp} - (gp'AVp)] / V, - ■ _^_JXTZ=^'- .,,^8 

f«.MC \^ 
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[029] It is also disclosed that all the above loadflow calculating systems 1) use the slack-start as 
the initial guess solution for increased processing acceleration, 2) skip from calculation, all the 
rows and columns corresponding to PV-nodes in factorized gain matrix [YV] while solving for 
[AV] in Q-V sub-problem, leading to some more processing acceleration, and 3) use the same 
indexing and addressing information for forming and storing both the constant gain matrices 
[Y9], and [YV] factorized using the same ordering regardless of node types leading to about 35% 
saving in computer storage/memory or alternatively calculating volume requirements. The 
SSDL-YY model is the simplest, easiest to implement and overall best in performance in terms of 
efficiency of calculations, reliability of convergence, and least computer memory requirements 
among all described in the following, and their simple variants with almost similar performance. 



[030 
[031 



[032 



(033 



[034 



[035 



BRIEF DESCRIPTION OF DRAWINGS 
Fig. 1 is a flow-chart of prior art FSDL method 

Fig. 2 is a flow-chart embodiment of the invented Super Super Decoupled Loadflow 
(SSDL-YY) calculation method 

Fig. 3 is a flow-chart embodiment of the invented Super Super Decoupled Loadflow 

(SSDL- BGX', SSDL-BGY, and SSDL-BGX) calculation methods 

Fig. 4 is a flow-chart embodiment of the invented Super Super Decoupled Loadflow 

(SSDL- X'GpvX') calculation method 

Fig. 5 is a flow-chart of the overall controlling method for an electrical power system 
involving loadflow calculation as a step which can be executed using one of the loadflow 
calculation methods embodied in Figs. 1, 2, 3 or 4 

Fig. 6 is a flow-chart of the simple special case of voltage control system in overall 
controlling system of Fig. 5 for an electrical power system 



v|036y Fjg/T ^s a one-line diagram of an exemplary 6-node powei^network h aving a 
\ rtj slack/sw ing/reference node, two PV-nodes, and three PQ-n odes 



DESCRIPTION OF A PREFERED EMBODYMENT 
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[037] A loadflow calculation is involved as a step in power flow control and/or voltage control 
in accordance with Fig. 5 or Fig. 6. A preferred embodiment of the present invention is described 
with reference to Fig. 7 as directed to achieving voltage control. ^^f(A . v^ . 
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[038] Fig. 7 is a simplified one-line diagram of an exemplary utility power network to which the 
present invention may be applied. The fundamentals of one-line diagrams are described in section 
6.11 of the text ELEMENTS OF POWER SYSTEM ANALYSIS, forth edition, by WiUiam D. 
Stevenson, Jr., McGrow-Hill Company, 1982. In Fig. 7, each thick vertical line is a network node. 
The nodes are interconnected in a desired manner by transmission lines and transformers each 
having its impedance, which appears in the loadflow models. Two transformers in Fig.7 are 
equipped with tap changers to control their turns ratios in order to control terminal voltage of 
node-1 and node-2 where large loads are connected. 

[039] Node-6 is a reference-node alternatively referred to as the slack or swing -node, 
representing the biggest power plant in a power network. Nodes-4 and -5 are PV-nodes where 
generators are connected, and nodes- 1, -2, and -3 are PQ-nodes where loads are connected. It 
should be noted that the nodes-4, -5, and -6 each represents a power plant that contains many 
generators in parallel operation. The single generator symbol at each of the nodes-4, -5, and -6 is 
equivalent of all generators in each plant. The power network further includes controllable circuit 
breakers located at each end of the transmission lines and transformers, and depicted by cross 
markings in one-line diagram of Fig. 7. The circuit breakers can be operated or in other words 
opened or closed manually by the power system operator or relevant circuit breakers operate 
automatically consequent of unhealthy or faulty operating conditions. The operation of one or 
more circuit breakers modify the configuration of the network. The arrows extending certain 
nodes represent loads. 

[040] A goal of the present invention is to provide a reliable and computationally efficient 
loadflow calculation that appears as a step in power flow control and/or voltage control systems 
of Fig. 5 and Fig. 6. However, the preferred embodiment of loadflow calculation as a step in 
control of terminal node voltages of PV-node generators and tap-changing transformers is 
illustrated in the flow diagram of Fig. 6 in which present invention resides in function steps 60 
and 62. 

[041] Short description of other possible embodiment of the present invention is also provided 
herein. The present invention relates to control of utility/industrial power networks of the types 
including plurality of power plants/generators and one or more motors/loads, and connected to 
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Other external utility. In the utility/industrial systems of this type, it is the usual practice to adjust 
the real and reactive power produced by each generator and each of the other sources including 
synchronous condensers and capacitor/inductor banks, in order to optimize the real and reactive 
power generation assignments of the system. Healthy or secure operation of the network can be 
shifted to optimized operation through corrective control produced by optimization functions 
without violation of security constraints. This is referred to as security constrained optimization 
of operation. Such an optimization is described in the United States Patent Number: 5,081,591 
dated Jan. 13, 1992: "Optimizing Reactive Power Distribution in an Industrial Power Network", 
where the present invention can be embodied by replacing the step nos. 56 and 66 each by a step 
of constant gain matrices [YG] and [YV], and replacing steps of "Exercise Newton-Raphson 
Algorithm" by steps of "Exercise Super Super Decoupled Loadflow Calculation" in places of 
steps 58 and 68. This is just to indicate the possible embodiment of the present invention in 
optimization functions like in many others including state estimation function. However, 
invention is being claimed through a simplified embodiment without optimization function as in 
Fig. 6 in this application. The inventive steps-60 and -62 in Fig.6 are different than those 
corresponding steps-56, and -58, which constitute a well known Newton-Raphson loadflow 
method, and were not inventive even in United States Patent Number: 5,081,591. 

[042] In Fig. 6, function step 10 provides stored impedance values of each network component 
in the system. This data is modified in a function step 12, which contains stored information 
about the open or close status of each circuit breaker. For each breaker that is open, the function 
step 12 assigns very high impedance to the associated line or transformer. The resulting data is 
than employed in a function step 14 to establish an admittance matrix for the power network. The 
data provided by function step 10 can be input by the computer operator from calculations based 
on measured values of impedance of each line and transformer, or on the basis of impedance 
measurements after the power network has been assembled. 

[043] Each of the transformers Tl and T2 in Fig. 7 is a tap changing transformer having a 
plurality of tap positions each representing a given transformation ratio. An indication of initially 
assigned transformation ratio for each transformer is provided by function step 20. 



« 
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[044] The indications provided by function steps 14, and 20 are supplied to a function step 60 in 
which constant gain matrices [Y0] and [YV] of any of the invented super decoupled loadflow 
mpdels are constructed, factorized and stored. The gain matrices [YG] and [YV] are conventional 
tools employed for solving Super Decoupled Loadflow model defined by equations (1) and (2) 
for a power system. 

[045] Indications of initial reactive power, or Q on each node, based on initial calculations or 
measurements, are provided by a function step 30 and these indications are used in function step 
32, to assign a Q level to each generator and motor. Initially, the Q assigned to each machine can 
be the same as the indicated Q value for the node to which that machine is connected. 

[046] An indication of measured real power, P, on each node is supplied by function step 40. 
Indications of assigned/specified/scheduled/set generating plant loads that are constituted by 
known program are provided by function step 42, which assigns the real power, P, load for each 
generating plant on the basis of the total P which must be generated within the power system. The 
value of P assigned to each power plant represents an economic optimum, and these values 
represent fixed constraints on the variations, which can be made by the system according to the 
present invention. The indications provided by function steps 40 and 42 are supplied to function 
step 44 which adjusts the P distribution on the various plant nodes accordingly. Function step 50 
assigns initial approximate or guess solution to begin iterative method of loadflow calculation, 
and reads data file of operating limits on power network components, such as maximum and 
minimum reactive power generation capability limits of PV-nodes generators. 

[047] The indications provided by function steps 32, 44, 50 and 60 are supplied to function step 
62 where inventive Fast Super Decoupled Loadflow calculation or Novel Fast Super Decoupled 
Loadflow calculation is carried out, the resuUs of which appear in function step 64. The loadflow 
calculation yields voltage magnitudes and voltage angles at PQ-nodes, real and reactive power 
generation by the slack/swing/reference node generator, voltage angles and reactive power 
generation indications at PV-nodes, and transformer turns ratio or tap position indications for tap 
changing transformers. The system stores in step 62 a representation of the reactive capability 
characteristic of each PV-node generator and these characteristics act as constraints on the 
reactive power that can be calculated for each PV-node generator for indication in step 64. The 
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indications provided in stop 64 actuate machine excitation control and transformer tap position 
control. All the loadflow calculation methods using SSDL models can be used to effect efficient 
and reliable voltage control in power systems as in the process flow diagram of Fig. 6. 

[048] Particular inventive loadflow calculation steps of all the SSDL methods are described 
followed by inventive SSDL models in terms of equations for determining elements of vectors 
[RP], [RQ], and elements of gain matrices [YG], and [YV] of equations (1) and (2), further 
followed by detailed steps of inventive loadflow calculation methods. Also, calculation 
accelerating and storage reducing techniques are described, which are not inventive. 

[049] The presence of transformed values of known/given/specified/scheduled/set quantities in 
the diagonal elements of the gain matrix [YV] of the Q-V sub-problem in equation (22), which 
takes different form for different manifestation of the SSDL-X'X' method, is brought about by 
such formulation of loadflow equations. The said transformed quantities in the diagonal 
elements of the gain matrix [YV] improved convergence and the reliability of obtaining 
converged loadflow calculation, when modified real power mismatch at any PQ-node-p is 
calculated as RPp = [APp' + (GppV Bpp') AQp'j / Vp^ which takes different form for different 
manifestation of the SSDL-X'X' method. 

[050] Rotation or transformation angle restriction to maximum of -48° in case of invented 
loadflow models for use in equations (35) and (36) or (27) and (28) is arrived at experimentally. 
The factor 0.9 in equations (38), and (39) is also determined experimentally. These said values of 
rotation angle -48° and the factor 0.9 are determined experimentally such that best possible 
convergence of loadflow calculation is obtained under wide range of network parameters and 
operating condition. 

[051] An altemate inventive system of loadflow calculation for Electrical Power System 
consists of, a Super Super Decoupled Loadflow: BGX'-version (SSDL-BGX') method, using 
simultaneous (IV, 19) iteration scheme, characterized in that 1) the modified real power 
mismatch at any PV-node-p is calculated as; 
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m 



RPp= [(APpA^p) - z GpqA Vq] / Vp _L nu) <n 



q=l 



and 2) the modified real power mismatch at any PQ-node-p is calculated as: 



m 



RPp = [(APpA^p) - 1 GpqAVq - (gp'AVp)] / Vp or , {Is) 



m 

RPp = [(APpA^p) - 1 GpqA Vq] / Vp 

q=l 



[052] A further altemate inventive system of loadflow calculation for Electrical Power System 
consists of, a Super Super Decoupled Loadflow:X'GpvX'-version (SSDL-X'GpyX') method, using 
simultaneous (IV, 10) iteration scheme, characterized in that 1) the modified real power 
mismatch at any PV-node-p is calculated as: 



m 




RPp = [(APpA^p) - S Gpq AVq] / (Kp* Vp) ^ - 

q=l 

and, 2) the modified real power mismatch at any PQ-node-p is calculated as: 
RPp = [{[APp' + (GppVBppOAQp'] / Vp} - (gp'AVp)] / Vp 

[053] The slack-start is to use the same voltage magnitude and angle as those of the 
slack/swing/reference node as the initial guess solution estimate for initiating the iterative 
loadflow calculation. With the specified/scheduled/set voltage magnitudes, PV-node voltage 
magnitudes are adjusted to their known values after the first P-9 iteration. This slack-start saves 
almost all effort of mismatch calculation in the first P-9 iteration. It requires only shunt flows 
fi-om each node to ground to be calculated at each node, because no flows occurs fi-om one node 
to another because they are at the same voltage magnitude and angle. 

[054] The same indexing and addressing information can be used for the storage of both the gain 
matrices [Y9] and [YV] as they are of the same dimension and sparsity structure as explained 
herein. This is achieved for all the inventive SSDL models. Voltage magnitude is 
specified/scheduled/set for PV-nodes where generators are connected, which each 
maintains/controls specified voltage magnitude at its node by changing its reactive power 
generation. Therefore, corrections in voltage magnitudes at PV-nodes are not required to be 
calculated as they do not change, and Q-V subproblem of equation (2) need to be solved only for 
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PQ-nodes voltage corrections. However, matrix [YV] is formed of dimension (m+k) x (m+k) 
including PV-nodes, and large diagonal value say, 10.0^^ is used for diagonal elements 
corresponding to PV-nodes to mask out the effect of off diagonal terms in the factorization 
process that normalizes the elements of a row in the matrix by its diagonal term. This process 
makes rows of PV-nodes in the gain matrix [YV] numerically absent despite physical presence. 
Because gain matrices [Y0] and [YV] are of the same power network, they are of the same 
dimension and sparsity structure. Therefore, they can be stored using the same indexing and 
addressing information leading to about 35% saving in^ computer memory or calculating volume 
requirements. Detailed description of the factorization process and storage schemes for sparse 
matrix is given by K.Zollenkopf, "Bi-factorization - Basic Calculational Algorithm and 
Programming Techniques", Large Sparse Sets of Linear Equations: Proceedings of Oxford 
conference of the Institute of Mathematics and its Application held in April 1 970, edited by 
J.K.Reid. 

[055] The efficiency of all inventive SSDL methods is increased by skipping all PV-nodes and 
factor elements with indices corresponding to PV-nodes while solving (5) for [AV]. In other 
words efficiency can be realized by skipping operations on rows/columns corresponding to PV- 
nodes in the forward-backward solution of (5). This is possible because rows and columns 
corresponding to PV-nodes are made numerically inactive by putting large value in diagonal 
terms of [YV]. This has been implemented and time saving of about 4% of the total solution time, 
including the time of input/output, could be realized. 

[056] Elements of gain matrices [Y9], and [YV] are possible to be defined independent of 
rotation or transformation angle in the inventive SSDL-YY model as explained herein. In super 
decoupling or transformation approach, rotation operators are applied to the complex node 
injections and the corresponding admittance values that relate them to the system state variables, 
transform the network equations such that branch or alternatively transmission line or transformer 
admittance appear to be almost entirely reactive. Better decoupling is thus realized. When rotated 
the complex branch admittance Ypq = Gpq + jBpq gets transformed into Ypq' = Gpq' + jBpq' , 
which is almost entirely reactive meaning Gpq' is almost of zero value. That means almost | Ypq' | 
= I Bpq' I . Therefore, by using | Ypq' | with the same algebraic sign as of Bpq' applied, in place of 
Bpq' in the definition of elements of gain matrices given by equations (38), (39), (59), (67), and 
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(82), the gain matrices [Y9], and [YV] are defined independent of rotation angles. This simplifies 
SSDL-YY model such that it is easy to program and implement, do not require performing steps 
involved in equations (102) to (108), and achieve improved acceleration in the loadflow 
calculation. 

Super Super Decoupled Loadflow: X'X'-version (SSDL-X'X') 

[057] The general model, in successive (10, IV) iteration scheme represented by sequence of 
equations (3) to (6), can be realized as SSDL-X'X', from which manifested are many versions. 
The elements of [RP], [RQ], [YQ] and [YV] are defined by (17) to (29). 

RPp = [APp' + (Gpp 7 Bpp') AQp'] / Vp^ -for PQ-nodes ( 1 7) 

RQp = [AQp' - (GppVBpp') APp'] / Vp -for PQ-nodes (18) 

RPp= [APp / (Kp*Vp2)] -forPV-nodes (19) 

YOpq = -1/Xpq' and YVp^ = -1/Xpq' (20) 

Yepp = E-Yepq and YVpp = bp' + I-YVp, (21) 

q>p q>p 

Where, bp' = -2bpCos<I)p or 
bp' = -bpCosOp + [QSHp' - (Gpp'/ Bpp') PSHp'] / Vs^ or 
bp' = 2[QSHp' - (Gpp'/ Bpp') PSHp'] / V^' (22) 

APp' = APpCosOp + AQpSin(I>p -for PQ-nodes (23) 

AQp' = AQpCosOp - APpSinOp -for PQ-nodes (24) 

PSHp' = PSHpCosOp + QSHpSinOp -for PQ-nodes (25) 

QSHp' = QSHpCosOp - PSHpSinOp -for PQ-nodes (26) 



CosOp = Absolute [Bpp / v/(Gpp^ + Bpp^)] > Cos (any angle from 0 to -90 degrees) (27) 
SinOp = -Absolute [Gpp / \/ (Gpp^* + B^] > Sin (any angle from 0 to -90 degrees) (28) 



Kp = Absolute (Bpp/YOpp) 



(29) 
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[058] Where, Gpp' and Bpp' are the transformed values of Gpp and Bpp as given by equation (1 06) 
in the appendix. The factor Kp of (29) is initially restricted to the minimum of 0.75 determined 
experimentally; however its restriction is lowered to the minimum value of 0.6 when its average 
over all less than 1.0 values at PV-nodes is less than 0.6. This factor is system and method 
independent. However it can be tuned for the best possible convergence for any given system. 
This statement is valid when the factor Kp is applied in the manner of equation (19) in all the 
models derived in the following from the most general model SSDL-X'X'. 

[059] The definition of YOpq in (20) is simplified because it does not explicitly state that it 
always takes the value of -Bpq for a branch connected between two PV-nodes or a PV-node and 
the slack-node. This fact should be understood implied in all the definitions of YGpq in this 
document. Xpq' in equation (20) is thejh-ans formed branch react ance defined in the appendix by 
equation (10 7).. 

[060] However, a whole new class of methods, corresponding to all those derived in the 
following and prior art, results when the factor Kp is used as a multiplier in the definition of RPp 
at PQ-nodes as in (30) instead of divider in RPp at PV-nodes as given in (19). This will cause 
changes only in equations (17), (19), and (20) as given in (30), (31), and (32). 

RPp = {[APp' + (Gpp' / Bpp') AQp'] / Vp^ } * Kp -for PQ-nodes (30) 
RPp= APp/Vp^ -for PV-nodes (31) 

Yepq=-Bpq and YVpq = -l/Xpq' (32) 

[061] Xpq' in equation (32) is the transformed branch reactance defined in the appendix by 
equation (107). The best performance of models of this new class has been realized when the 
factor Kp, applied in a manner of equation (30) leading to changes as in (30) to (32), is 
unrestricted. That means it can take any value as given by equation (29). 

Super Super Decoupled Loadflow: YY-version (SSDL-YY) 

[062] If unrestricted rotation is applied to complex branch admittance and transformed 
susceptance is taken as admittance magnitude value with the same algebraic sign and transformed 
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conductance is assumed zero, the SSDL-X'X' method reduces to SSDL-YY. Though, this model 
is not very sensitive to the restriction applied to nodal rotation angles, SSDL-YY presented here 
restricts rotation angles to the maximum of -48 degrees determined experimentally for the best 
possible convergence from non linearity considerations, which is depicted by equations (35) and 
(36). However, it gives closely similar performance over wide range of restriction applied to the 
nodal rotation angles say, from -36 to -90 degrees. , 



RPp = (APpCos<Dp + AQpSin<Dp)/ Vp^ -forPQ-nodes T-y, HuT (33) 

RQp = (AQpCosOp - APpSinOp )A^p -for PQ-nodes (34) 

Cos(I>p = Absolute (Bpp/^/(Gpp^ + Bpp^) > Cos (-48°) '.^"^r^tf^S^^ L. /..^^^'^ 

SinOp = -Absolute (Gpp / /(G^T+B^) > Sin (-48°) 0^ (im\^ C3 h) ^^r) ^ 

cl4/cii pb'''vi , — 

RPp = APp / (KpVp^) -forPV-nodes niMM^e, i"^'^ (37) 

-Ypq -for branch r/x ratio < 3.0 'T\^-m\y ^. Tt, HmV 

- (Bp, + 0.9(Ypc-Bp<,)) -for branch r/x ratio > 3.0 ^^i) 1- ^ '^^^'^'^ ^""^ ' 

-Bpq -for branches connected between two PV-nodes or a 

PV-node and the slack-node (38) 



YGnn — 



YVp, 



-Ypq -for branch r/x ratio < 3.0 

- (Bpq + 0.9(Ypq-Bpq)) -for branch r/x ratio > 3.0 (39) 



Y0pp = I-YGpq and YVpp = bp' + Z-YVpq (40) 

q>p q>p 

bp* = (QSHpCosOp - PSHpSinOp / Vs^) - bpCosOp or 



bp'=2(QSHpCosOp - PSHpSinOp )/ Vs^ ^Jcp^^iW J (41) 




[063] where, Kp is defined in equation (29), which is initially restricted to the minimum 
value of 0.75 determined experimentally; however its restriction is lowered to the minimum 
value of 0.6 when its average over all less than 1.0 values at PV nodes is less than 0.6. 
Restrictions to the factor Kp as stated in the above is system independent. However it can be tuned 
for the best possible convergence for any given system. In case of systems of only PQ-nodes 
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and without any PV-nodes, equations (38) and (39) simply be taken as YOpq = YVpq = -Ypq. 
Branch admittance magnitude in (38) and (39) is of the same algebraic sign as its susceptance. 
Elements of the two gain matrices differ in that diagonal elements of [YV] additionally contain 
the b' values given by equations (40) and (41) and in respect of elements corresponding to 
branches connected between two PV-nodes or a PV-node and the slack-node. Equations (35) and 
(36) with inequality sign implies that nodal rotation angles are restricted to maximum of -48 
degrees for SSDL-YY. The model consists of equations (3) to (6), (33) to (41), and (29). In two 
simple variations of the SSDL-YY model, one is to make YYr,n=Y%n and the other is to make 



|064l) The steps of loadflow calculation method, SSDL-YY method are shown in the flowchart 



similar letters in the following. Double lettered steps are the characteristic steps of SSDL-YY 
method. The words "Read system data" in Step-a correspond to step- 10 and step-20 in Fig. 5, and 
step- 14, step-20, step-32, step-44, step-50 in Fig. 6. AH other steps in the following correspond to 
step-30 in Fig.5, and step-60, step-62, and step-64 in Fig. 6. 

a. Read system data and assign an initial approximate solution. If better solution estimate is 
not available, set voltage magnitude and angle of all nodes equal to those of the slack- 
node. This is referred to as the slack-start. 

b. Form nodal admittance matrix, and Initialize iteration count ITRP = ITRQ= r = 0 

cc. Compute Cosine and Sine of nodal rotation angles using equations (35) and (36), and store 
them. If they, respectively, are less than the Cosine and Sine of -48 degrees, equate them, 
respectively, to those of -48 degrees. 

dd. Form (m+k) x (m+k) size matrices [YG] and [YV] of (1) and (2) respectively each in a 
compact storage exploiting sparsity. The matrices are formed using equations (38), (39), 
(40), and (41). In [YV] matrix, replace diagonal elements corresponding to PV-nodes by 
very large value (say, 10.0*^). In case [YV] is of dimension (m x m), this is not required to 
be performed. Factorize [Y0] and [YV] using the same ordering of nodes regardless of 
node-types and store them using the same indexing and addressing information. In case 
[YV] is of dimension (m x m), it is factorized using different ordering than that of [Y0]. 





of Fig. 2. Referring to the flowchart of Fig.2, different steps are elaborated in steps marked with 
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e. Compute residues [AP] at PQ- and PV-nodes and [AQ] at only PQ-nodes. If all are less 

than the tolerance (e), proceed to step-n. Otherwise follow the next step, 
ff. Compute the vector of modified residues [RP] as in (33) for PQ-nodes, and using (37) and 
(29) for PV-nodes. 

g. Solve (3) for [AG] and update voltage angles using, [9] = [9] + [A9]. 

h. Set voltage magnitudes of PV-nodes equal to the specified values, and Increment the 
iteration count ITRP=ITRP+1 and i=(ITRP+ITRQ)/2. 

i. Compute residues [AP] at PQ- and PV-nodes and [AQ] at PQ-nodes only. If all are less 
than the tolerance (e), proceed to step-n. Otherwise follow the next step. 

j. Compute the vector of modified residues [RQ] as in (34) for only PQ-nodes. 
k. Solve (5) for [AV] and update PQ-node magnitudes using [V] = [V] + [AV]. While 
solving equation (5), skip all the rows and columns corresponding to PV-nodes. 
Calculate reactive power generation at PV-nodes and tap positions of t£^-changing 
transformers. If the maximum and minimum reactive power generation capability and 
transformer tap position limits are violated, implement the violated physical limits and 
adjust the loadflow solution by the method described in "LTC Transformers and MVAR 
violations in the Fast Decoupled Load Flow, IEEE Trans., PAS-101, No.9, PP. 3328- 
3332, September 1982". 

J, ^ m. Increment the iteration count ITRQ=ITRQ+1 and r=(ITRP+ITRQ)/2, and Proceed to 

step-e. 




7r> —f^- From calculated values of voltage magnitude and voltage angle at PQ-nodes, voltage 

angle and reactive power generation at PV-nodes, and tap position of tap changing 
transformers, calculate power flows through power network components. 

Super Super Decoupled Loadflow: XX-verslon (SSDL-XX) 

[065] If no or zero rotation is apphed, the SSDL-X'X' method reduces to SSDL-XX, which is 
the simplest form of SSDL-X'X'. The SSDL-XX model comprises equations (3) to (6), (42) to 
(48), and (29). 

RPp = [APp + (Gpp / Bpp) AQp] / Vp2 -for PQ-nodes (42) 

RQp = [AQp - (Gpp / Bpp) APp] / Vp -for PQ-nodes (43) 
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RPp = APp / (KpVp^) -for PV-nodes (44) 



Y6pq = - 1 .0/Xpq -for all other branches 

-Bpq -for branches connected between two PV-nodes or a 

PV-node and the slack-node (45) 
YVpq= -1.0/Xpq -for all branches (46) 

YGpp = Z- YGpq and YVpp = bp' + Z-YVpq (47) 

q>p q>p 



bp' = -2bp or 
bp' = -bp + [QSHp - (Gpp / Bpp) PSHp]A^s^ or 

bp' = 2[QSHp - (Gpp/ Bpp) PSHplA^s"" (48) 

[066] where, Kp is defined in equation (29). This is the simplest method with very good 
performance for distribution networks in absence of PV-nodes or alternatively for systems 
containing only PQ-nodes. The large value of the difference [(1/X)-B], particularly for high R/X 
ratio branches connected to PV-nodes, creates modeling error when PV-nodes are present in a 
system. 

Super Super Decoupled Loadflow: BX-version (SSDL-BX) 

[067] If super decoupling is applied only to QV-sub problem, the SSDL-XX model reduces to 
SSDL-BX, which makes it perform better for systems containing PV-nodes. The SSDL-BX 
model comprises equations (3) to (6), (49) to (51), (47) and (48). This model can be referred to as 
Advanced BX-Fast Decoupled Loadflow. 

RPp= APp/Vp^ -for all nodes (49) 

RQp = [AQp - (Gpp / Bpp) APp] / Vp -for PQ-nodes (50) 

YBpq = -Bpq and Y Vpq = - 1 /Xpq (5 1 ) 



[068] It should be noted that Amerongen's General-purpose Fast Decoupled Loadflow model 

described in "A general-purpose version of the Fast Decoupled Loadflow", IEEE Transactions, 

PWRS-4, pp.760-770, May 1989; has turned out to be an approximation of this model. The 
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approximation involved is only in equation (50). However, numerical performance is found to be 
only slightly better but more reliable than that of the Amerongen's method. 

Super Super Decoupled Loadflow: X'B'-version (SSDL-X'B') 

[069] This model consists of equations (3) to (6), (52) to (57), and (23) to (29). 
RPp=[APp' + (Gpp7Bpp')AQp']/Vp2 -forPQ-nodes (52) 

RQp= AQpVVp -forPQ-nodes (53) 

RPp= [APp / (Kp* Vp')] -forPV-nodes (54) 

Yep<,= -l/Xpq' and YVpq = -Bpq' (55) 

Y9pp = Z-Yepc and YVpp = bp' + I-YVp, (56) 

Where, bp' = -2bpCos<Dp or 

bp' = -bpCos<Dp + QSHpV Vs^ or 

bp' = 2QSHpV Vs' (57) 

[070] Where, APp', AQp', PSHp', QSHp', CosOp, SinOp, Kp are defined in (23) to (29), and Bpq' 
and Xpq' are the transformed branch susceptance and reactance defined in the appendix by 
equations (105) and (107). Best performance of this model could be achieved by restricting Op in 
(27) and (28) to less than or equal to -48°. 

Super Super Decoupled Loadflow: YB'-version (SSDL-YB') 

(071 J The equation (52) in SSDL-X'B' implies unrestricted % is applied and it can take values 
up to -90 degrees. Therefore, (52) can be modified to (58) with consequent modification of (55) 
into (59). This model consists of equations (3) to (6), (58), (53), (54), (59), (56) and (57), and (23) 
to (29). 



2 

P 



RPp = [APp* Absolute. [Bpp / ^(Gpp^ + Bpp^)] + AQp* [-Absolute [Gpp / /(Gpp^ + Bpp^ / V, 

-for PQ-nodes (58) 
YOpq = -Ypq and YVp^ = -Bpq' (59) 
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[072] The best performance of this method could be achieved by restricting C>p in (27) and (28) 
to less than or equal to -48 degrees. Where, APp' , AQp' , PSHp' , QSHp' , Cos<Dp , S'm% , Kp are 
defined in (23) to (29), and Bpq' is the transformed branch susceptance defined in the appendix by 
equations ( 1 05). 

Super Super Decoupled Loadflow: B'X'-version (SSDL-B'X') 

[073] This model consists of equations (3) to (6), (60) to (65), and (23) to (29). 



RPp = APpV Vp^ -for PQ-nodes (60) 

RQp = [AQp'-(GppVBpp') APp']/Vp -for PQ-nodes (61) 

RPp= [APp / (Kp*Vp^)] -forPV-nodes (62) 

Yepq=-Bpq' and YVpq = -l/Xpq' (63) 

Yepp = E-Yepq and YVpp = bp' + Z-YVpq (64) 

q>p q>p 

Where, bp' = -2bpCos(I)p or 

bp' = -bpCosOp + [QSHp' - (GppV Bpp') PSHp']A^s^ or 

bp' = 2[QSHp' - (GppV Bpp') PSHp'lA^s' (65) 



[074] Where, APp', AQp', PSHp', QSHp', CosOp, SinOp, Kp are defined in (23) to (29), and Bpq' 
and Xpq' are the transformed branch susceptance and reactance defined in the appendix by 
equations (105) and (107). Best performance of this model could be achieved by restricting Op in 
(27) and (28) to less than, equal to -48°. 

Super Super Decoupled Loadflow: B'Y-version (SSDL- B'V) 

[075] The realtion (61) in SSDL-B'X' implies unrestricted Op is applied and it can take values 

up to -90 degrees. Therefore, (61) can be modified to (66) with consequent modification of (63) 
into (67). 

RQp = [AQp* Absolute [Bpp / \/(Gpp^ + Bpp^)] - APp* [-Absolute [Gpp / v^G^7+^^]] / Vp^ 

-for PQ-nodes (66) 

YOpq = -Bpq' and YVpq = -Ypq (67) 
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[076] This model consists of equations (3) to (6), (60), (66), (62), (67), (64) and (65), and (23) to 
(29). Best performance of this model could be achieved by restricting <Dp in (27) and (28) to less 
than or equal to -48 degrees. Where, APp' , AQp' , PSHp' , QSHp' , CosOp , Sin<Dp , Kp are defined 
in (23) to (29), and Bpq' is the transformed branch susceptance defined in the appendix by 
equation (105). 

Simultaneous (IV, 10) Iteration Scheme 

[077] An ideal to be approached for the decoupled loadflow models is the constant matrix 
loadflow model described in "Fast Super Decoupled Loadflow", lEE proceedings Part-C, 
Vol.139, No.l, pp. 13-20, January 1992; referred in this document as BGGB-model. In an attempt 
to imitate it, a decoupled class of models with simultaneous (IV, 19) iteration scheme depicted by 
sequence of equations (68) to (72) is invented. This scheme involves only one mismatch 
calculation in an iteration. The correction vector is calculated in two separate parts of voltage 
magnitude and voltage angle corrections without intermediate updating. Each iteration involves 
one calculation of [RQ], [AV], and [RP], [AG] to update [V] and [9]. 

[AV] = [YV]-'[RQ] (68) 
[RP] = [APA^] - [G] [AV] (69) 
[A9] = [Y9] ■' [RP] (70) 
[9] = [9] + [A9] (71) 
[V] = [V] + [AV] (72) 

[078] In this invented class, each model differs only in the definition of elements of [RQ] and 
[YV]. The accuracy of models depends only on the accuracy of calculation of [AV]. The greater 
the angular spread of branches terminating at PQ-nodes, the greater the inaccuracy in the 
calculation of [AV]. 

Super Super Decoupled Loadflow: BGX'-version (SSDL-BGX') 

[079] Numerical performance could further be improved by organizing the solution in a 
simultaneous (IV, 19) iteration scheme represented by sequence of equations (68) to (72). The 
elements of [RP], [RQ], [Y0] and [YV] are defined by (73) to (80). 



RQp = [AQp' - (Gpp' / Bpp') APp'] / Vp 
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A 



m 



RPp= [(APp/Vp).ZGpqAV,]/Vp 



-for all nodes 



or 



(74) 



m 



YGpq — B 



YGpp = Y-YQ 



q>p 



RPp = [(APpA^p) - 1 GpqAVq - (gp'AVp)] / Vp -for PQ-nodes, and 

-forPV-nodes f^P^^H76) 
and YVpq = -l/Xpq' m)j jg£^j^^^ (77) 

and YVoB = bo' + Z-YVoa^'''*^ T'^M^'c^e*/ '^Jf ^(yj^^^^'^i^ 



RPp = [(APpA^p) - 1 GpqAVq] / Vp 

q=l 



pq 



pq 



pp v/p ■ Tpq 
q>p 



Where, 



bp' = -ZbpCosOp 



bp' = - bpCosOp + [QSHp' - (GppVBpp') PSHp'JA^s^ 
bp' = 2[QSHp' - (Gpp'/Bpp') PSHp']/ Vs^ 



or 
or 



(79) 




gp' = 2bpSin<Dp or 
gp' = bpSinOp + [PSHp' + (GppV Bpp')QSHp'] A^s^ or ^ 

gp' = 2[PSHp' + (GppV Bpp')QSHp']A^s^ 

. § ^- [080] Where, APp' , AQp' , PSHp' , QSHp' , Cos<Dp , SinOp are defined in (23) to (28), andXpq] 
"X. I is the transf orme d branch reactance defined in the appendix by equations (107) . The SSDL-BGX' 
model comprises equations (68) to (80), and (23) to (28); Best possible convergence could be 
achieved by restricting rotation angles Op in the range of -lO" to -20° in equations (27) and (28). 
The model is embodied in the flow-chart of Fig.3. 

Super Super Decoupled Loadflow: BGY-version (SSDL-BGY) 



[081] If unrestricted rotation is applied and transformed susceptance is taken as admittance 
values and transformed conductance is assumed zero, the SSDL-BGX' model reduces to SSDL- 
BGY as. defined by equations (81), (74), (75), (76), (82) to (85). ; &ci yiX<i^ 



RQp = AQp'A^p = (AQpGosOp - APpSin(Dp)A^p 



-for PQ-nodes 



Y0pq -Bpq 



and 



YV„n = -Y„„ 
1 V pq I pq 
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Yepp = Z-Yepq and YVpp = bp' + Z-YVpq (83) 

Where, bp' = -2bpCos4)p or 



7 




bp' = -bpCos<Dp + (QSHpCosOp - PSHpSin<Dp)A's^ or ^ j-^ 

bp' = 2(QSHpCos<Dp - PSHpSin(Dp)A^s^ (84) ^ 

gp' = 2bpSin(Dp or 
gp' = bpSinOp + (PSHpCosOp + QSHpSin<Dp)A^s^ or 
bp' = 2(PSHpCosOp + QSHpSina)p)AVs^ 

[082] The SSDL-BGY model comprises equations (68) to (72), and (81), (74), (75), (76), (82) to 
(85). It is the special case of the SSDL-BGX' model. 

Super Super Decoupled Loadflow: BGX-version (SSDL-BGX) 

[083] If no or zero rotation is applied, the SSDL-BGX' model reduces to SSDL-BGX as defined 
by equations (86), (74), (75), (76), (87) to (90). ^ ^\ 

RQp = [AQp - (Gpp / Bpp)APp] / Vp -for PQ-nodes -otwt i, ii| (86) r J^y^h^ 

YGpq = -Bpq and YVp^ = -l/Xp, to/>0' C<t\) y-h h'^e^f «PJ/ "^-n ^^^j 

YGpp = I-Y0pq and YVpp = bp' + |-YVpq (88) /<' «i tr) 

bp' = -2bp or = -bp+[QSHp -(Gpp/Bpp)PSHp]A^s' or = 2[QSHp-(Gpp/Bpp)PSHp]/Vs' (S9j)L^h\-r) 
gp' = 0.0 or = [PSHp +(Gpp/Bpp)QSHp]A^s' or = 2[PSHp +(Gpp/Bpp)QSHp]A^s' (m) ajd dj 

[084] The SSDL-BGX model comprises equations (68) to (72), and (86), (74), (75), (76), (87) to 
(90)_It is again the special case of the SSDL-BGX' model. 

[OSSTiThe stens of Inadflow nalr.nlfltinn mpthnHc <5<snT _Rr;Y' ccr^T nr;v ocrir un.v^''^^<l^'f^A~~ 
methods are shown in the flowchart of Fig. 3. Referring to the flowchart of Fig.3, different steps 
are elaborated in steps marked with similar letters in the following. Triple lettered steps are the 
characteristic steps of SSDL-BGX', SSDL-BGY and SSDL-BGX methods. The words "Read 
system data" in Step-a correspond to step- 10 and step-20 in Fig. 5, and step- 14, step-20, step-32, 
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step-44, step-50 in Fig. 6. All other steps in the following correspond to step-30 in Fig.5, and 
step-60, step~62, and step-64 in Fig. 6. 

a. Read system data and assign an initial approximate solution. If better solution estimate is 
not available, set voltage magnitude and angle of all nodes equal to those of the slack- 
node. This is referred to as the slack-start 

b. Form nodal admittance matrix, and Initialize iteration count ITR = 0, 

ccc. Compute Sine and Cosine of nodal rotation angles using equations (28) and (27), and store 
them. If they, respectively, are less than the Sine and Cosine of any angle set say, in the 
range -10 to -20 degrees, equate them, respectively, to those of the same angle in the range 
-10 to -20 degrees. In case of zero rotation, Sine and Cosine value vectors are not required 
to be stored, 

ddd. Form (m+k) x (m+k) size matrices [Y9] and [YV] of (1) and (2) respectively each in a 
compact storage exploiting sparsity 

1) In case of SSDL-BGX'-method, the matrices are formed using equations (77), (78), 
and (79) 

2) In case of SSDL-BGY-method, the matrices are formed using equations (82), (83), 
and (84) 

3) In case of SSDL-BGX-method, the matrices are formed using equations (87), (88), 
and (89) 

In [YV] matrix, replace diagonal elements corresponding to PV-nodes by very large value 
say, 10.0'". In case [YV] is of dimension (m x m), this is not required to be performed. 
Factorize [Y9] and [YV] using the same ordering of nodes regardless of node-types and 
store them using the same indexing and addressing information. In case [YV] is of 
dimension (m x m), it is factorized using different ordering than that of [YO]. 

e. Compute residues [AP] at PQ- and PV-nodes and [AQ] at only PQ-nodes. If all are less 
than the tolerance (e), proceed to step-n. Otherwise follow the next step. 

fff. Compute the vector of modified residues [RQ] using (73) in case of SSDL-BGX' , using 
(81) in case of SSDL-BGY, and using (86) in case of SSDL-BGX method for only PQ- 
nodes. Solve (68) for [AV]. While solving equation (68), skip all the rows and columns 
corresponding to PV-nodes. Compute the vector of modified residues [RP] using (74) or 
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ggg- 




m. 



n. 



(75) and (76). Solve (70) for [AG]. 

Update voltage angles using, [0] = [0] + [AG], and update PQ-node voltage magnitudes 
using [V] = [V] + [AV]. Set voltage magnitudes of PV-nodes equal to the specified values. 
Calculate reactive power generation at PV-nodes and tap positions of tap-changing 
transformers. If the maximum and minimum reactive power generation capability and 
transformer tap position limits are violated, implement the violated physical limits and 
adjust the loadflow solution by the method described in "LTC Transformers and MVAR 
violations in the Fast Decoupled Load Flow, IEEE Trans., PAS-101, No.9, PP. 3328-3332, 
September 1982". 

Increment the iteration count ITRQ=ITRQ+1 and r=(ITRP+ITRQ)/2, and Proceed to 
step-e. 

From calculated values of voltage magnitude and voltage angle at PQ-nodes, voltage angle 
and reactive power generation at PV-nodes, and tap position of tap changing transformers, 
calculate power flows through power network components. 



1086] The SSDL-BGX', SSDL-BGY and SSDL-BGX methods differ only in steps-ccc and -ddd 
defining gain matrices, and step-fff for calculating [RP] and [RQ]. 

Super Super Decoupled Loadflow: X'GpvX'-version (SSDL- X'GpyX') 

[087] Numerical performance could be improved by organizing the solution in a simultaneous 
(IV, 10) iteration scheme represented by sequence of equations (68) to (72). The elements of 
[RP], [RQ], [YG] and [YV] for this model are defined by (91) to (94), and (78), (79), and (80). 



RQp = [AQp' - (Gpp' /Bpp')APp'] / V, 



m 



RPp= [(APpA^p)-EGpqAVq]/(Kp*Vp) 

q=l 



YGpq= -1/Xpq' 



and 



YVpq = -l/Xpq' 



-for PQ-nodes 



RPp = [{[APp' + (Gpp7Bpp')AQp'] / Vp} - (gp'AVp)] / Vp -for PQ-nodes 
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[088) Where, APp' , AQp' , PSHp' , QSHp' , Cos% , SinOp , Kp are defined in (23) to (29), and 
Xpq' is the transformed branch reactance defined in the appendix by the equation (107). Again, if 
unrestricted rotation is applied and transformed susceptance is taken as admittance values and 
transformed conductance is assumed zero, the SSDL-X'GpvX' model reduces to SSDL-YGpvY. 
If no or zero rotation is applied, the SSDL-X'GpyX' model reduces to SSDL-XGpyX. The SSDL- 
X'GpvX' model comprises equations (68) to (72), (91) to (94), (78) to (80), and (23) to (29). 

[0891 The steps of loadflow calculation methods, SSDL-X'GpvX', SSDL-YGpvY and SSDL- 
XGpvX methods are shown in the flowchart of Fig. 4. Referring to the flowchart of Fig.4, 
different steps are elaborated in steps marked with similar letters in the following. Four lettered 
steps are the characteristic steps of SSDL-X'GpyX', SSDL-YGpyY and SSDL-XGpyX methods. 
These methods are useful particularly for distribution systems without PV-nodes. The words 
"Read system data" in Step-a correspond to step- 10 and step-20 in Fig. 5, and step-14, step-20, 
step-32, step-44, step-50 in Fig. 6. All other steps in the following correspond to step-30 in Fig.5, 
and step-60, step-62, and step-64 in Fig. 6. 

a. Read system data and assign an initial approximate solution. If better solution estimate is 
not available, set voltage magnitude and angle of all nodes equal to those of the slack- 
node. This is referred to as the slack-start. 

b. Form nodal admittance matrix, and Initialize iteration count ITR= 0. 

cccc. Compute Sine and Cosine of nodal rotation angles using equations (28) and (27), store 
them. If they, respectively, are less than the Sine and Cosine of any angle set say, 0 to -90 
degrees, equate them, respectively, to those of the same angle in the range 0 to -90 
degrees. In case of zero rotation, Sine and Cosine vectors are not required to be stored. 

dddd. Form (m+k) x (m+k) size matrices [Y6] and [YV] of (1) and (2) respectively each in a 
compact storage exploiting sparsity using equations (94), (78), and (79). 
In [YV] matrix, replace diagonal elements corresponding to PV-nodes by very large value 
say, 10.0^^. In case [YV] is of dimension (m x m), this is not required to be performed. 
Factorize [Y9] and [YV] using the same ordering of nodes regardless of node-types and 
store them using the same indexing and addressing information. In case [YV] is of 
dimension (m x m), it is factorized using different ordering than that of [YO]. 

e. Compute residues [AP] at PQ- and PV-nodes, and [AQ] at only PQ-nodes. If all are less 
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than the tolerance (e), proceed to step-n. Otherwise follow the next step, 
ffff. Compute [RQ] using (9 1 ) for only PQ-nodes. Solve (68) for [AV]. While solving equation 
(68), skip all the rows and columns corresponding to PV-nodes. Compute the vector of 
modified residues [RP] using equations (92), (93), (80) and (29). Solve (70) for [AG], 
ggg. Update voltage angles using, [0] = [9] + [AG], and update PQ-node voltage magnitudes 

using [V] = [V] + [AV]. Set volta ge magnitudes of PV-nodes equal to the specified values. fJ-^eJ 




Calculate reactive power generation at PV-nodes and tap positions of tap-changing 
(xM^ transformers. If the maximum and minimum reactive power generation capability and 

transformer tap position limits are violated, implement the violated physical limits and 
adjust the loadflow solution by the method described in "LTC Transformers and MVAR 
violations in the Fast Decoupled Load Flow, IEEE Trans., PAS- 101, No.9, PP. 3328-3332, 
September 1982". 

I ^m. Increment the iteration count ITRQ=ITRQ+1 and r=(ITRP+ITRQ)/2, and Proceed to 

step-e. 

-jt) ^ n. From calculated values of voltage magnitude and voltage angle at PQ-nodes, voltage angle 

and reactive power generation at PV-nodes, and tap position of tap changing transformers, 
calculate power flows through power network components. 

MODE^USING GLOBAL CORRECTIONS 

[090] It is possible that all the above described models can be organized to produce corrections 
to the initial estimate solution. It involves storage of the vector of modified residues and replacing 
the equations (17), (18), (19), which are the most general form of equations appearing in the 
description, and (49) by (95), (96), (97), and (98) respectively, and (4) or (71) and (6) or (72) 
respectively by (99) and (100). Superscript '0' in equations (99) and (100) indicates the initial 
solution estimate. y ./ £^ cj^^ ^ ^ 

RPp^ = [(APp^)' + (GppV Bpp') (AQpT]/(Vpf + RPp('-') (95) i<{i^\^) 

RQp^- [(AQpO' - (GppVBpp') (APpT]/(VpO' + RQp^^ '^ (96) ^ 

RP/ = APp' / [Kp(Vpf ] + RP/-') (97) i^iM 1 ) 

RPp^= APpV(VpO^ + RPp(^-') 
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v;=Vp°+AVp' 



C^ht»)_^(99) 



RECTANGULAR COORDINATE FORMULATIONS OF LOADFLOW MODELS 

[091] This involves following changes in the equations describing the loadflow models 
formulated in polar coordinates. 

(i) Replace G and AG respectively by f and Af in equations (1), (3), (4), (70), (71) and (99) 

(ii) Replace V and AV respectively by e and Ae in equations (2), (5), (6), (68), (69), (72) and 
(100) 

(iii) Replace Vp by Cp and Vs by Cs in equations (17) to (19), (22), (30), (31), (33), (34), (37), 
(41) to (44), (48) to (50), (52) to (54), (57), (58), (60) to (62), (65), (66), (73), (74), (75), 
(76), (79) to (81), (84) to (86), (89) to (93), (95), (96), (97), (98). The subscript ^s' indicates 
the slack-node variable. 

(iv) After calculation of corrections to the imaginary part of complex voltage (Af) of P V-nodes 
and updating the imaginary component (f) of PV-nodes, calculate real component by: 

ep =-NiVsp'-fp' f^-rfk)^m) 



Vp interchangeably in all expressions of RPp, and the use of 1.0 for in all expressions of bp' 
and gp'. However, in equations (74), (75), (76), (92), and (93) interchangeable use of Vpand Vp^ is 
allowed only when any one of the dividing term Vp assumes value of 1 .0. 

[093] While the description above refers to particular embodiments of the present invention, it 
will be understood that many modifications may be made without departing from the spirit 
thereof The accompanying claims are intended to cover such modifications as would fall within 
the true scope and spirit of the present invention. 

[094] The presently disclosed embodiments are therefore to be considered in all respect as 
illustrative and not restrictive, the scope of the invention being indicated by the appended claims 
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in addition to the foregoing description, and all changes which come within the meaning and 
range of equivalency of the claims are therefore intended to be embraced therein. 



APPENDIX 



[095] The following steps give the procedure for determining transformed/rotated branch 
admittance necessary in forming symmetrical gain matrices [YV] and/or [Y9] of the SSDL 
models. 

1 . Compute; Op = arctan (Gpp/Bpp) and 

<Dq = arctan (GVBW f^j, rA)^(102) 

2. Compute the average of rotations at the terminal nodes p and q of a branch: 

3. Compare Oav with the Limiting Rotation Angle (LRA) and let Oav to be the smaller of the 

two: 

Oav = minimum (Oav, LRA) Qc^VYk) 04) 

4. Compute transformed pq-th element of the admittance matrix: 

Ypq' = Gpq' + jBpq' = (CoS Oav + jSin Oav) (Gpq + jBpq) (105) 

Vpp-Gpp'.jBpp'=E-(Gp,'.jBp,') 

5. Note that the transformed branch reactance is: 

Xpq' = Bpq V(Gpq'2 + Bpq'^) Q ^(1 07) 

Similarly, Xpp' = Bpp7(Gpp'^ + Bpp'^) {\0^tir) ^^^^^S) 
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CLAIMS 



The embodiments of the invention in which an exclusive property or privilege is claimed are 
defined as follows: 



A method of forming/defining and solving a model of the power network to effect control 
of voltages and power flows in a power system, comprising the steps of: 




obtaining on-line/simulated data of open/close status of switches and circuit breakers in 
the power network, and reading data of operating limits of the network 
components including PV-node, a generator-node where Real-Power-P and 
Voltage-Magnitude-V are given/assigned/specified/set, generators maximum and 
minimum reactive power generation capability limits and transformers tap position 
lim its. 



obtaining on-line readings of given/assigned/specified/set real and reactive power at PQ- 
nodes, the load-nodes where Real-Power-P and Reactive-Power-Q are 
given/assigned/specified/set, real power and voltage magnitude at PV-nodes, 
voltage magnitude and angle at the reference/slack node, and transformer turns 
ratios, which are the controlled variables/parameters. 



initiating loadflow calculation with initial approximate/guess solution of the same voltage 
magnitude and angle as those of the slack/reference node for all the other nodes 

referred to as the slack-start, 



i forming and storing factorized gain matrices [Y0] and [YV] using the same indexing and 

addressing information for both as they are of the same dimension and sparsity 

structure, wherein said [Y0] relate vector of modified real power mismatches [RP] 

^ to angle corrections vector [AG] in equation [RP] = [YG] [AG] referred to as P-0 

sub-problem, and said [YV] relate vector of modified reactive power mismatches 

[RQ] to voltage magnitude corrections vector [AV] in equation [RQ] = [YV] [AV] 

referred to as Q-V sub-problem, 
I " ■ — 

^J> /restricting transformation/rotation angle % to maximum -48'' in determining transformed 
j^^V^ V ^ ^^^^ reactive power mismatch as, 

. 0> APp'= APpCosOp 4- AQpSinOp -for PQ-nodes (23) 

^ - A 

AQp'= AQpCosOp - APpSinOp -for PQ-nodes (24) 

Wherein, APp and AQp are real and reactive power mismatches at node-p. 




N% 'v'^ 
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calculating modified real and reactive power mismatches as given in the following in the 
most general form of equations that take different form for different Super Super 
Decoupled Loadflow model: 



RPp = [APp' + (GppVBpp') AQp'] A^p2 
RQp = [AQp' - (GppVBpp') APp'] / Vp 



-for PQ-nodes 
-for PQ-nodes 



(17) 
(18) 



and calculating modified real power mismatch at a PV-node as, 



-for PV-nodes 



(19) 
(29) 



RPp= APp/(KpVp^ ) 

Wherein, Kp = Absolute (BppA^Gpp) 
and Vp is voltage magnitude at node-p, and Bpp is imaginary part of the diagonal 
element Ypp of the admittance matrix without network shunts, and YOpp is the 
diagonal element of the gain matrix [Y9], 




using network shunt parameter bp' that appears in diagonal elements of gain matrix [YV] 
as given in the following in the most general form of equations that take different 
form for different loadflow model: 



bp' = -bpCosOp+ [QSHp' - (GppVBpp') PSHp'] / Vs^ 
bp' = 2[QSHp' - (Gpp'/Bpp') PSHp'] / V^' 



or 



(22) 



wherein, Gpp'and Bpp' are the real and imaginary parts of the transformed diagonal 
element Ypp' of the admittance matrix without network shunts, bp is network shunt 



at node-p, Vj is slack-node voltage magnitude, and 



PSHp' = PSHpCosOp + QSHpSinOp 
QSHp' = QSHpCosa)p- PSHpSinOp 



-for PQ-nodes 
-for PQ-nodes 



(25) 
(26) 




A 



wherein, PSHp, and QSHp are given/specified/scheduled/set real and reactive 
power respectively, 

performing loadflow calculation by solving a super decoupled loadflow model of the 
power network defined by set of equations [RP] = [YG] [AG] and [RQ] = [YV] 
[AV] employing successive (19, IV) iteration scheme, wherein each iteration 
involves one calculation of [RP] and [AG] to update voltage angle vector [G] and 
then one calculation of [RQ] and [AV] to update voltage magnitude vector [V], to 
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A 



.V 



A 



calculate values of the voltage angle and the voltage magnitude at PQ-nodes,. 
voltage angle and reactive power generation at PV-nodes, and turns ratio of tap- 
changing transformers in dependence on the set of said obtained-online readings, 
or given/scheduled/specified/set values of controlled variables/parameters and 
physical limits of operation of the network components, 
evaluating loadflow calculation for any of the over loaded power network components and 

for under/over voltage at any of the network nodes, 
Correcting one or more controlled parameters and repeating the calculating, performing, 
evaluating, and correcting steps until evaluating step finds no over loaded 
^components and no under/over voltages in the power network, and 
Nfeffecting a change in the power flowing through network components and voltage 
0^ magnitudes and angles at the nodes of the power network by actually 

implementing the finally obtained values of controlled variables/parameters after 
evaluating step finds a good power system or alternatively a power network 
without any overloaded components and under/over voltages, which finally 
obtained controlled variables/parameters however are stored in case of simulation 
for acting upon fast in case the simulated, event actually occurs. 





2. A 



method as defined in claim- 1 wherein loadflow calculation involving formation and 
solution of super super decoupled loadflow model, employing simultaneous (IV, 
10) iteration scheme is characterized in that it involve only one time calculation of 
real and reactive power mismatches in an iteration along with modified real power 
mismatch calculation, depending on super super decoupled loadflow model used, 



either by: 



m 



RPp = [(APpA^p) . I Gpq AV J / Vp 



-for all nodes 



or 



(74) 



m 



RPp = [(APpAVp) - Z GpqAVq - (gp'AVp)] / Vp -for PQ-nodes, and 

q=l 



(75) 



m 



\ 



RPp - [(APp/Vp) - Z GpqAVq] / Vp 

q=l 

OR 



-for PV-nodes 



(76) 
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RPp = [ {[APp' + (GppV Bpp')AQp'] / Vp} - (gp'AVp)] / Vp -for PQ-nodes (92) 



m 



RPp= [(APpA^p) - 1 Gp, AV,] / (Kp*Vp) -forPV-nodes (93) 



q=l 



wherein, m is the number of PQ-nodes, and the symborA' preceding any variable 
represents mismatch or correction in the variable. 

A simple system for controlling generator and transformer voltages in an electrical power 
utiUty containing plurality of electromechanical rotating machines, transformers 
and electrical loads connected in a network, each machine having a reactive power 
characteristic and an excitation element which is controllable for adjusting the 
reactive power generated or absorbed by the machine, and some of the 
transformers each having a tap changing element which is controllable for 
adjusting turns ratio or alternatively terminal voltage of the transformer, said 
system comprising: 

means defining and solving one of the Super Super Decoupled Loadflow models of the 
power network characterized in claim- 1 or claim-2 for providing an indication of 
the quantity of reactive power to be supplied by each generator including the 
slack/reference node generator, and for providing an indication of turns ratio of 
each tap-changing transformer in dependence on the set of obtained-online 
readings or given/specified/set controlled network variables/parameters, and 
physical Umits of operation of the network components, 
machine control means connected to the said means defining and solving one of the Super 
Super Decoupled Loadflow models of the power network and to the excitation 
elements of the rotating machines for controlling the operation of the excitation 
elements of machines to produce or absorb the amount of reactive power indicated 
by said means defining and solving one of the Super Super Decoupled Loadflow 
models of the power network in dependence on the set of obtained-online readings 
or given/specified/set controlled network variables/parameters, and physical limits 
of excitation elements, 

transformer tap position control means connected to the said means defining and solving 
one of the Super Super Decoupled Loadflow models of the power network and to 
the tap changing elements of the controllable transformers for controlling the 
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operation of the tap changing elements to adjust the turns ratios of transformers 
indicated by the said means defining and solving Super Super Decoupled 
Loadflow model of the power network in dependence on the set of obtained-online 
readings or given/specified/set controlled network variables/parameters, and 
operating limits of the tap-changing elements. 
A system as defined in claim-3 wherein the power network includes a plurality of nodes 
each connected to at least one of: a slack/reference generator; a rotating machine; 
and an electrical load, and the said means defining and solving one of the Super 
Super Decoup;ed Loadflow models of the power network receives representations 
of selected values of the real and reactive power flow firom each machine and to 
each load, and the model is operative for producing calculated values for the 
reactive power quantity to be produced or absorbed by each machine. 
A system as defined in claim-4 wherein the power network further has at least one 
transformer having an adjustable transformer turns ratio, and the said means 
defining and solving one of the Super Super Decoupled Loadflow models of the 
power network is further operative for producing a calculated value of the 
transformer transformation/turns ratio. 
A system as defined in claim-3 wherein said machine control means are connected to said 
excitation element of each machine for controlling the operation of the excitation 
element of each machine, and wherein said transformer turns ratio control means 
are connected to said transformer tap changing element of each transformer for 
controlling the operation of the tap changing element of each transformer. 
A method for controlling generator and transformer voltages in an electrical power utility 
containing plurality of electromechanical rotating machines, transformers and 
electrical loads connected in a network, each machine having a reactive power 
characteristic and excitation element which is controllable for adjusting the 
reactive power generated or absorbed by the machine, and some of the 
transformers each having tap changing element which is controllable for adjusting 
turns ratio or alternatively terminal voltage of the transformer, said method 
comprising: 
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creating and solving any of the said Super Super Decoupled Loadflow models of the 
power network as characterized in claim- 1 or claim-2 for providing an indication 
of turns ratio of each of the tap-changing transformers and the quantity of reactive 
power to be supplied by each generator in dependence on the set of obtained- 
online readings or given/specified/set controlled network variables/parameters, and 
physical limits of operation of the network components, - 

controlling the operation of the excitation elements of machines to produce or absorb the 
amount of reactive power, and controlling tap changing elements of transformers 
to adjust transformer turns ratio indicated by means creating/forming/defming and 

7 

solving any of the said Super Super Decoupled Loadflow models as characterized 
in claim- 1 or claim-2 in dependence on the set of obtained-onliiie readings or 
given/specified/set controlled network variables/parameters<( and physical limits of 
operation of the network components. 

^ Q ^ 8. A method as defined in claim-7 wherein said step of controlling is carried out-to control 

the excitation element of each machine. 

^ \ ^ 9. A method as defined in claim-7 wherein said step of controlling is carried out to control 

the tap-changing element of each controllable transformer. 
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SYSTEM OF SUPER SUPER DECOUPLED LOADFLOW CALCULATION FOR 

ELECTRICAL POWER SYSTEM 

ABSTRACT 

A method of performing loadflow calculations for controlling voltages and power flow in a 
power network by reading on-line data of given/specified/scheduled/set network 
variables/parameters and using control means, so that no component of the power network is 
overloaded as well as there is no over/under voltage at any nodes in the network following a small 
or large disturbances. The invented generahzed Super Super Decoupled Loadflow (SSDL) 
calculation method is characterized in that 1) modified real power mismatch at any PQ-node-p is 
calculated as RPp = [APp' + (GppV Bpp') AQp'] / Vp^ which takes different form for different 
manifestation of the generalized version SSDL-X'X' method, 2) transformed ^values of 
known/given/specified/scheduled/set quantities in the diagonal elements of the gain matrix [YV] 
of the Q- V sub-problem are present, and 3) transformation angles are restricted to maximum of - 
48° particularly for the most successful version SSDL-YY method, and these inventive loadflow 
calculation steps together yield some processing acceleration and consequent efficiency gains, 
and are each individually inventive. The other two Super Super Decoupled Loadflow methods: 
BGX* version (SSDL-BGX') and X^GpvX' version (SSDL-X'GpvX^ are characterized in the use of 
simultaneous (IV, 16) iteration scheme thereby calculating mismatches only once in each 
iteration and consequent efficiency gain. 
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Fig. 7: An Exemplary 6-node Power System 

Nodes: 1, 2, 3 are PQ-nodes 

Nodes: 4 and 5 are PV-nodes 

Nodes: 6 is the slack/swing/reference node 

Transformers Tl and T2 are tap-changing 



